Metabolic priming by multiple enzyme systems supports glycolysis, HIF1α stabilisation, and human cancer cell survival in early hypoxia.

Adaptation to chronic hypoxia occurs through changes in protein expression, which are controlled by hypoxia-inducible factor 1α (HIF1α) and are necessary for cancer cell survival. However, the mechanisms that enable cancer cells to adapt in early hypoxia, before the HIF1α-mediated transcription programme is fully established, remain poorly understood. Here we show in human breast cancer cells, that within 3 h of hypoxia exposure, glycolytic flux increases in a HIF1α-independent manner but is limited by NAD+ availability. Glycolytic ATP maintenance and cell survival in early hypoxia rely on reserve lactate dehydrogenase A capacity as well as the activity of glutamate-oxoglutarate transaminase 1 (GOT1), an enzyme that fuels malate dehydrogenase 1 (MDH1)-derived NAD+. In addition, GOT1 maintains low α-ketoglutarate levels, thereby limiting prolyl hydroxylase activity to promote HIF1α stabilisation in early hypoxia and enable robust HIF1α target gene expression in later hypoxia. Our findings reveal that, in normoxia, multiple enzyme systems maintain cells in a primed state ready to support increased glycolysis and HIF1α stabilisation upon oxygen limitation, until other adaptive processes that require more time are fully established.

A branching model of lineage differentiation underpinning the neurogenic potential of enteric glia.

Glial cells have been proposed as a source of neural progenitors, but the mechanisms underpinning the neurogenic potential of adult glia are not known. Using single cell transcriptomic profiling, we show that enteric glial cells represent a cell state attained by autonomic neural crest cells as they transition along a linear differentiation trajectory that allows them to retain neurogenic potential while acquiring mature glial functions. Key neurogenic loci in early enteric nervous system progenitors remain in open chromatin configuration in mature enteric glia, thus facilitating neuronal differentiation under appropriate conditions. Molecular profiling and gene targeting of enteric glial cells in a cell culture model of enteric neurogenesis and a gut injury model demonstrate that neuronal differentiation of glia is driven by transcriptional programs employed in vivo by early progenitors. Our work provides mechanistic insight into the regulatory landscape underpinning the development of intestinal neural circuits and generates a platform for advancing glial cells as therapeutic agents for the treatment of neural deficits.

Functional cross-talk between allosteric effects of activating and inhibiting ligands underlies PKM2 regulation.

Several enzymes can simultaneously interact with multiple intracellular metabolites, however, how the allosteric effects of distinct ligands are integrated to coordinately control enzymatic activity remains poorly understood. We addressed this question using, as a model system, the glycolytic enzyme pyruvate kinase M2 (PKM2). We show that the PKM2 activator fructose 1,6-bisphosphate (FBP) alone promotes tetramerisation and increases PKM2 activity, but addition of the inhibitor L-phenylalanine (Phe) prevents maximal activation of FBP-bound PKM2 tetramers. We developed a method, AlloHubMat, that uses eigenvalue decomposition of mutual information derived from molecular dynamics trajectories to identify residues that mediate FBP-induced allostery. Experimental mutagenesis of these residues identified PKM2 variants in which activation by FBP remains intact but cannot be attenuated by Phe. Our findings reveal residues involved in FBP-induced allostery that enable the integration of allosteric input from Phe and provide a paradigm for the coordinate regulation of enzymatic activity by simultaneous allosteric inputs.

The TAB1-p38α complex aggravates myocardial injury and can be targeted by small molecules.

Inhibiting MAPK14 (p38α) diminishes cardiac damage in myocardial ischemia. During myocardial ischemia, p38α interacts with TAB1, a scaffold protein, which promotes p38α autoactivation; active p38α (pp38α) then transphosphorylates TAB1. Previously, we solved the X-ray structure of the p38α-TAB1 (residues 384-412) complex. Here, we further characterize the interaction by solving the structure of the pp38α-TAB1 (residues 1-438) complex in the active state. Based on this information, we created a global knock-in (KI) mouse with substitution of 4 residues on TAB1 that we show are required for docking onto p38α. Whereas ablating p38α or TAB1 resulted in early embryonal lethality, the TAB1-KI mice were viable and had no appreciable alteration in their lymphocyte repertoire or myocardial transcriptional profile; nonetheless, following in vivo regional myocardial ischemia, infarction volume was significantly reduced and the transphosphorylation of TAB1 was disabled. Unexpectedly, the activation of myocardial p38α during ischemia was only mildly attenuated in TAB1-KI hearts. We also identified a group of fragments able to disrupt the interaction between p38α and TAB1. We conclude that the interaction between the 2 proteins can be targeted with small molecules. The data reveal that it is possible to selectively inhibit signaling downstream of p38α to attenuate ischemic injury.

Erratum: Genome editing reveals a role for OCT4 in human embryogenesis.

This corrects the article DOI: 10.1038/nature24033.

Genome editing reveals a role for OCT4 in human embryogenesis.

Despite their fundamental biological and clinical importance, the molecular mechanisms that regulate the first cell fate decisions in the human embryo are not well understood. Here we use CRISPR-Cas9-mediated genome editing to investigate the function of the pluripotency transcription factor OCT4 during human embryogenesis. We identified an efficient OCT4-targeting guide RNA using an inducible human embryonic stem cell-based system and microinjection of mouse zygotes. Using these refined methods, we efficiently and specifically targeted the gene encoding OCT4 (POU5F1) in diploid human zygotes and found that blastocyst development was compromised. Transcriptomics analysis revealed that, in POU5F1-null cells, gene expression was downregulated not only for extra-embryonic trophectoderm genes, such as CDX2, but also for regulators of the pluripotent epiblast, including NANOG. By contrast, Pou5f1-null mouse embryos maintained the expression of orthologous genes, and blastocyst development was established, but maintenance was compromised. We conclude that CRISPR-Cas9-mediated genome editing is a powerful method for investigating gene function in the context of human development.

An engineered photoswitchable mammalian pyruvate kinase.

Changes in allosteric regulation of glycolytic enzymes have been linked to metabolic reprogramming involved in cancer. Remarkably, allosteric mechanisms control enzyme function at significantly shorter time-scales compared to the long-term effects of metabolic reprogramming on cell proliferation. It remains unclear if and how the speed and reversibility afforded by rapid allosteric control of metabolic enzymes is important for cell proliferation. Tools that allow specific, dynamic modulation of enzymatic activities in mammalian cells would help address this question. Towards this goal, we have used molecular dynamics simulations to guide the design of mPKM2 internal light/oxygen/voltage-sensitive domain 2 (LOV2) fusion at position D24 (PiL[D24]), an engineered pyruvate kinase M2 (PKM2) variant that harbours an insertion of the light-sensing LOV2 domain from Avena Sativa within a region implicated in allosteric regulation by fructose 1,6-bisphosphate (FBP). The LOV2 photoreaction is preserved in the PiL[D24] chimera and causes secondary structure changes that are associated with a 30% decrease in the Km of the enzyme for phosphoenolpyruvate resulting in increased pyruvate kinase activity after light exposure. Importantly, this change in activity is reversible upon light withdrawal. Expression of PiL[D24] in cells leads to light-induced increase in labelling of pyruvate from glucose. PiL[D24] therefore could provide a means to modulate cellular glucose metabolism in a remote manner and paves the way for studying the importance of rapid allosteric phenomena in the regulation of metabolism and enzyme control.

Lineage-dependent spatial and functional organization of the mammalian enteric nervous system.

The enteric nervous system (ENS) is essential for digestive function and gut homeostasis. Here we show that the amorphous neuroglia networks of the mouse ENS are composed of overlapping clonal units founded by postmigratory neural crest-derived progenitors. The spatial configuration of ENS clones depends on proliferation-driven local interactions of ENS progenitors with lineally unrelated neuroectodermal cells, the ordered colonization of the serosa-mucosa axis by clonal descendants, and gut expansion. Single-cell transcriptomics and mutagenesis analysis delineated dynamic molecular states of ENS progenitors and identified RET as a regulator of neurogenic commitment. Clonally related enteric neurons exhibit synchronous activity in response to network stimulation. Thus, lineage relationships underpin the organization of the peripheral nervous system.

A Gene Regulatory Network Balances Neural and Mesoderm Specification during Vertebrate Trunk Development.

Transcriptional networks, regulated by extracellular signals, control cell fate decisions and determine the size and composition of developing tissues. One example is the network controlling bipotent neuromesodermal progenitors (NMPs) that fuel embryo elongation by generating spinal cord and trunk mesoderm tissue. Here, we use single-cell transcriptomics to identify the molecular signature of NMPs and reverse engineer the mechanism that regulates their differentiation. Together with genetic perturbations, this reveals a transcriptional network that integrates opposing retinoic acid (RA) and Wnt signals to determine the rate at which cells enter and exit the NMP state. RA, produced by newly generated mesodermal cells, provides feedback that initiates NMP generation and induces neural differentiation, thereby coordinating the production of neural and mesodermal tissue. Together, the data define a regulatory network architecture that balances the generation of different cell types from bipotential progenitors in order to facilitate orderly axis elongation.

Coevolved Mutations Reveal Distinct Architectures for Two Core Proteins in the Bacterial Flagellar Motor.

Switching of bacterial flagellar rotation is caused by large domain movements of the FliG protein triggered by binding of the signal protein CheY to FliM. FliG and FliM form adjacent multi-subunit arrays within the basal body C-ring. The movements alter the interaction of the FliG C-terminal (FliGC) "torque" helix with the stator complexes. Atomic models based on the Salmonella entrovar C-ring electron microscopy reconstruction have implications for switching, but lack consensus on the relative locations of the FliG armadillo (ARM) domains (amino-terminal (FliGN), middle (FliGM) and FliGC) as well as changes during chemotaxis. The generality of the Salmonella model is challenged by the variation in motor morphology and response between species. We studied coevolved residue mutations to determine the unifying elements of switch architecture. Residue interactions, measured by their coevolution, were formalized as a network, guided by structural data. Our measurements reveal a common design with dedicated switch and motor modules. The FliM middle domain (FliMM) has extensive connectivity most simply explained by conserved intra and inter-subunit contacts. In contrast, FliG has patchy, complex architecture. Conserved structural motifs form interacting nodes in the coevolution network that wire FliMM to the FliGC C-terminal, four-helix motor module (C3-6). FliG C3-6 coevolution is organized around the torque helix, differently from other ARM domains. The nodes form separated, surface-proximal patches that are targeted by deleterious mutations as in other allosteric systems. The dominant node is formed by the EHPQ motif at the FliMMFliGM contact interface and adjacent helix residues at a central location within FliGM. The node interacts with nodes in the N-terminal FliGc α-helix triad (ARM-C) and FliGN. ARM-C, separated from C3-6 by the MFVF motif, has poor intra-network connectivity consistent with its variable orientation revealed by structural data. ARM-C could be the convertor element that provides mechanistic and species diversity.

The Phylogenetic Signature Underlying ATP Synthase c-Ring Compliance.

The proton-driven ATP synthase (FOF1) is comprised of two rotary, stepping motors (FO and F1) coupled by an elastic power transmission. The elastic compliance resides in the rotor module that includes the membrane-embedded FO c-ring. Proton transport by FO is firmly coupled to the rotation of the c-ring relative to other FO subunits (ab2). It drives ATP synthesis. We used a computational method to investigate the contribution of the c-ring to the total elastic compliance. We performed principal component analysis of conformational ensembles built using distance constraints from the bovine mitochondrial c-ring x-ray structure. Angular rotary twist, the dominant ring motion, was estimated to show that the c-ring accounted in part for the measured compliance. Ring rotation was entrained to rotation of the external helix within each hairpin-shaped c-subunit in the ring. Ensembles of monomer and dimers extracted from complete c-rings showed that the coupling between collective ring and the individual subunit motions was independent of the size of the c-ring, which varies between organisms. Molecular determinants were identified by covariance analysis of residue coevolution and structural-alphabet-based local dynamics correlations. The residue coevolution gave a readout of subunit architecture. The dynamic couplings revealed that the hinge for both ring and subunit helix rotations was constructed from the proton-binding site and the adjacent glycine motif (IB-GGGG) in the midmembrane plane. IB-GGGG motifs were linked by long-range couplings across the ring, while intrasubunit couplings connected the motif to the conserved cytoplasmic loop and adjacent segments. The correlation with principal collective motions shows that the couplings underlie both ring rotary and bending motions. Noncontact couplings between IB-GGGG motifs matched the coevolution signal as well as contact couplings. The residue coevolution reflects the physiological importance of the dynamics that may link proton transfer to ring compliance.

The crystal structure of the human titin:obscurin complex reveals a conserved yet specific muscle M-band zipper module.

M10 is the most C-terminal immunoglobulin (Ig) domain of the giant protein titin and a frequent target of disease-linked mutations. Currently, it is the only known muscle Ig domain able to interact with two alternative ligands-obscurin and obscurin-like-1 (Obsl1)-in different sarcomeric subregions. Obscurin and Obsl1 use their homologous N-terminal Ig domain (O1 in obscurin and OL1 in Obsl1) to bind M10 in a mutually exclusive manner. We present here the X-ray structure of the human titin:obscurin M10:O1 complex extending our previous work on the M10:OL1 interaction. Similar to M10:OL1, the M10:O1 complex displays a chevron-shaped antiparallel Ig-Ig architecture held together by a conserved molecular interface, which we validated by isothermal titration calorimetry and sorting experiments in neonatal rat cardiomyocytes. O1, although structurally related to OL1 and M10, both members of the intermediate set (I-set) Ig family, presents an intriguing switch of its ßA' strand. This leads to structural differences between the complexes, particularly for the "open side" of the chevron-shaped assembly. A bioinformatics analysis reveals that the ßA'-switch observed for O1 is rare and that it is involved in mediating protein-protein interactions. Molecular dynamics simulations also suggest that this topological alteration substantially increases local flexibility compared to the conventional I-set Ig domains. The O1/OL1 Ig domains are candidate discriminatory structural modules potentially directing the binding of specific additional partners at the M-band. Cellular sorting experiments in neonatal rat cardiomyocytes are consistent with the view that the titin:obscurin/Obsl1 complexes might be a platform for higher-order interactions.

In vitro generation of neuromesodermal progenitors reveals distinct roles for wnt signalling in the specification of spinal cord and paraxial mesoderm identity.

Cells of the spinal cord and somites arise from shared, dual-fated precursors, located towards the posterior of the elongating embryo. Here we show that these neuromesodermal progenitors (NMPs) can readily be generated in vitro from mouse and human pluripotent stem cells by activating Wnt and Fgf signalling, timed to emulate in vivo development. Similar to NMPs in vivo, these cells co-express the neural factor Sox2 and the mesodermal factor Brachyury and differentiate into neural and paraxial mesoderm in vitro and in vivo. The neural cells produced by NMPs have spinal cord but not anterior neural identity and can differentiate into spinal cord motor neurons. This is consistent with the shared origin of spinal cord and somites and the distinct ontogeny of the anterior and posterior nervous system. Systematic analysis of the transcriptome during differentiation identifies the molecular correlates of each of the cell identities and the routes by which they are obtained. Moreover, we take advantage of the system to provide evidence that Brachyury represses neural differentiation and that signals from mesoderm are not necessary to induce the posterior identity of spinal cord cells. This indicates that the mesoderm inducing and posteriorising functions of Wnt signalling represent two molecularly separate activities. Together the data illustrate how reverse engineering normal developmental mechanisms allows the differentiation of specific cell types in vitro and the analysis of previous difficult to access aspects of embryo development.

Design and application of implicit solvent models in biomolecular simulations.

We review implicit solvent models and their parametrisation by introducing the concepts and recent devlopments of the most popular models with a focus on parametrisation via force matching. An overview of recent applications of the solvation energy term in protein dynamics, modelling, design and prediction is given to illustrate the usability and versatility of implicit solvation in reproducing the physical behaviour of biomolecular systems. Limitations of implicit modes are discussed through the example of more challenging systems like nucleic acids and membranes.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet.

MOTIVATION: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies. AVAILABILITY: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools.

Decrypting Prion Protein Conversion into a ß-Rich Conformer by Molecular Dynamics.

Prion diseases are fatal neurodegenerative diseases characterized by the formation of ß-rich oligomers and the accumulation of amyloid fibrillar deposits in the central nervous system. Understanding the conversion of the cellular prion protein into its ß-rich polymeric conformers is fundamental to tackling the early stages of the development of prion diseases. In this paper, we have identified unfolding and refolding steps critical to the conversion into a ß-rich conformer for different constructs of the ovine prion protein by molecular dynamics simulations. By combining our results with in vitro experiments, we show that the folded C-terminus of the ovine prion protein is able to recurrently undergo a drastic conformational change by displacement of the H1 helix, uncovering of the H2H3 domain, and formation of persistent ß-sheets between H2 and H3 residues. The observed ß-sheets refold toward the C-terminus exposing what we call a "bending region" comprising residues 204-214. This is strikingly coincident with the region harboring mutations determining the fate of the prion oligomerization process. The ß-rich intermediate is used here for the construction of a putative model for the assembly into an oligomeric aggregate. The results presented here confirm the importance of the H2H3 domain for prion oligomer formation and therefore its potential use as molecular target in the design of novel prion inhibitors.

Specific Conformational States of Ras GTPase upon Effector Binding.

To uncover the structural and dynamical determinants involved in the highly specific binding of Ras GTPase to its effectors, the conformational states of Ras in uncomplexed form and complexed to the downstream effectors Byr2, PI3Kγ, PLCε, and RalGDS were investigated using molecular dynamics and cross-comparison of the trajectories. The subtle changes in the dynamics and conformations of Ras upon effector binding require an analysis that targets local changes independent of global motions. Using a structural alphabet, a computational procedure is proposed to quantify local conformational changes. Positions detected by this approach were characterized as either specific for a particular effector, specific for an effector domain type, or as effector unspecific. A set of nine structurally connected residues (Ras residues 5-8, 32-35, 39-42, 55-59, 73-78, and 161-165), which link the effector binding site to the distant C-terminus, changed dynamics upon effector binding, indicating a potential effector-unspecific signaling route within the Ras structure. Additional conformational changes were detected along the N-terminus of the central ß-sheet. Besides the Ras residues at the effector interface (e.g., D33, E37, D38, and Y40), which adopt effector-specific local conformations, the binding signal propagates from the interface to distant hot-spot residues, in particular to Y5 and D57. The results of this study reveal possible conformational mechanisms for the stabilization of the active state of Ras upon downstream effector binding and for the structural determinants responsible for effector specificity.

Implicit Solvation Parameters Derived from Explicit Water Forces in Large-Scale Molecular Dynamics Simulations.

Implicit solvation is a mean force approach to model solvent forces acting on a solute molecule. It is frequently used in molecular simulations to reduce the computational cost of solvent treatment. In the first instance, the free energy of solvation and the associated solvent-solute forces can be approximated by a function of the solvent-accessible surface area (SASA) of the solute and differentiated by an atom-specific solvation parameter σ(i) (SASA). A procedure for the determination of values for the σ(i) (SASA) parameters through matching of explicit and implicit solvation forces is proposed. Using the results of Molecular Dynamics simulations of 188 topologically diverse protein structures in water and in implicit solvent, values for the σ(i) (SASA) parameters for atom types i of the standard amino acids in the GROMOS force field have been determined. A simplified representation based on groups of atom types σ(g) (SASA) was obtained via partitioning of the atom-type σ(i) (SASA) distributions by dynamic programming. Three groups of atom types with well separated parameter ranges were obtained, and their performance in implicit versus explicit simulations was assessed. The solvent forces are available at http://mathbio.nimr.mrc.ac.uk/wiki/Solvent_Forces.

Urea-Water Solvation Forces on Prion Structures.

Solvation forces are crucial determinants in the equilibrium between the folded and unfolded state of proteins. Particularly interesting are the solvent forces of denaturing solvent mixtures on folded and misfolded states of proteins involved in neurodegeneration. The C-terminal globular domain of the ovine prion protein (1UW3) and its analogue H2H3 in the α-rich and ß-rich conformation were used as model structures to study the solvation forces in 4 M aqueous urea using molecular dynamics. The model structures display very different secondary structures and solvent exposures. Most protein atoms favor interactions with urea over interactions with water. The force difference between protein-urea and protein-water interactions correlates with hydrophobicity; i.e., urea interacts preferentially with hydrophobic atoms, in agreement with results from solvent transfer experiments. Solvent Shannon entropy maps illustrate the mobility gradient of the urea-water mixture from the first solvation shell to the bulk. Single urea molecules replace water in the first solvation shell preferably at locations of relatively high solvent entropy.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics.

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.